IP GAS-HE-6/95 



^-matrix description of the finite-temperature 

noneqilibrium media 

J.ManjavidzeS 



Abstract 

The paper contains the real-time perturbation theory for description of a statistical 
system with the nonuniform temperature distribution. The formalism based on the Wigner- 
functions approach. The perturbation theory is formulated in terms of the local-temperature 
Green functions. 
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1 Introduction 



The aim of this article is to construct the perturbation theory for generating functional 
of Wigner functions |T], for the case of nonuniform temperature distribution. As an 
example of interesting system one can have in mind the process of a very large number of 
hadrons creation at the high-energy collisions. The phenomenology of high-multiplicity 
processes was given in jf]. [|. 

In terms of QCD it is the rear process of the cold quark-gluon plasma formation. Qual- 
itatively this is a process of total dissipation of high kinetic energy density initial state. 
One can consider this process also as the process of dissipation of a high-temperature 
local fluctuation in a low temperature equilibrium media. Generally speaking, the tem- 
perature freely evolves in this processes and is not distributed uniformly at least on the 
early stages. In this paper we will consider the general theory of such processes which 
can be used not only in the particles physics. 

We should adopt the standard 5"-matrix formalism which is applicable to any nonequi- 
librium processes. In this microcanonical approach the temperature T will be introduced 
as the Lagrange multiplier and the physical (measurable) value of T will be defined by 
the equation of states at the very end of calculations. Using standard terminology ||, 
we will deal with the "mechanical" perturbations only and it will not be necessary to 
divide the perturbations on "thermal" and "mechanical" ones || [| (see also Sec. 5). 

The usual Kubo-Martin-Schvinger (KMS) periodic boundary conditions |7], [H| can not 
be applied here since they are applicable for the equilibrium case only [M (see also [I2|). 
We will introduce the boundary conditions "by hands" , modeling the environment of the 
system [|l^] . Supposing that the system is in a vacuum we will have usual field-theoretical 
vacuum boundary condition (Sec. 2). We will consider also the system in the background 
field of black-body radiation (Sec. 4). Last one restores the theory with KMS boundary 
condition in the equilibrium limit. 

Calculating the generating functional of Wigner functions the local temperature dis- 
tribution will be introduced: T(x,t) = l//3(x,t) is the temperature in the measurement 
point (x,t) (Sec. 3). In other words, we will divide the "measuring device" (but not the 
system as usually was done, [113]) on the cells the dimension of which tends to zero. The 
differential measure Df3(x, t) will be defined taking into account the energy-momentum 
conservation law. 



2 Vacuum boundary condition 

The probability r(P) of in- into out-states transition with fixed total 4-momentum P can 
be calculated using the n- into m-particle transition amplitude a n ^ m : 

n,m 'Tl- m - J 

n m 

x5 (4) (P _ J- q k )5^(P -J2 Pk )\a n , m \ 2 i (2.1) 



where 

n n j3„ 

fe=i k=i y zn ) ze \ik) 



Eq. ( |2.1| ) is the basic formula of our calculations. The microcanonical description was 



introduced in [|T^] considering the Fourier transformation of 5-functions. 



The amplitude a n . m looks as follows [ |T2]j : 

n m 

a n , m ((q)n, (p)m) = Y[ 4>(qk) II }*(pk)Z((f>), (2.3) 

k=l k=l 

where qk{Pk) are the momentum of in(out)-going particles and the annihilation operator 

0(g) = J d 4 xe-^^x), = j^-, (2.4) 

was introduced. Correspondingly, <p*{p) is the creation operator. One can put the auxil- 
iary field (f>(x) equal to zero at the end of calculation. The vacuum into vacuum transition 
amplitude in presence of external field <fi 

Z(<p) = J D$ e «c + (*)-iV 0+ (#+^) ( 2 5 ) 

is defined on the Mills' complex time contour C + [|TJ[], i.e. C + : t — > t + ie,e > 0. In 
eq.( |2.5| ) Sc + is the free part of the action and Vc + describes the interactions. 
In this section we will propose the vacuum boundary condition: 

da^d^ = (2.6) 



where o"oo is the infinitely far hypersurface. 

We start consideration from the assumption that the temperature fluctuations are 
large scale. In a cell the dimension of which is much smaller then the fluctuation scale of 
temperature we can assume that the temperature is a "good" parameter. (The "good" 
parameter means that the corresponding fluctuations are Gaussian.) 

Let us surround the interaction region, i.e. the system under consideration, by N cells 
with known space-time position and let us propose that we can measure the energy and 
momentum of groups of in- and out-going particles in each cell. The 4-dimension of cells 
can not be arbitrary small in this case because of the quantum uncertainty principle. 

To describe this situation we decompose 5-functions in ( [2.1| ) on the product of (N+ 1) 
^-functions: 



n r N n v N 

6^(p-j2q k )= / nw*-E?M)}* (4, (^-EW. ( 2 - 7 ) 

k=l J v=\ k=l v=l 

where qu yU are the momentum of k-th in-going particle in the z/-th cell and Q u is the total 
4-momenta of n u in-going particles in this cell. The same decomposition will be used for 



the second (^-function in ( j2.ll ). Inserting this decompositions into ( j2.1| ) we must take into 
account the multinomial character of particles decomposition on N groups. This will give 
the coefficient: 



n\ „ , " , m! 



M n ~E n ^) i ' :5 K (m-J2 m »)i ( 2 - 8 ) 



ni!---njv! ,7^ mi\---m N \ 



where is the Kronecker's 5-function. 
In result, the quantity 



r((Q) N , (P) N ) = J2 J \a(n,m)\ 2 x 



(n.m) 

x mn ^^^ (4) (^ - e ^ n - ep v )} (2.9) 

t/=l fe=l n v k=l k=l k=l 

describes the probability that in the z/-th cell we measure the fluxes of in-going particles 
with total 4-momentum Q v and of out-going particles with the total 4- momentum P v . 
The sequence of this two measurements is not fixed. 

The Fourier transformation of (^-functions in ( ^9| ) gives the formula: 

r((Q) N ,(P) N ) = / II (2 ^r ^^ ^^-^^^CC"-)^, Kk), (2.10) 
where i?((a_)jv, (a+)jv) = -R(a_ 5 i, a-,2---, ct-,v; a+,2, ck+,aO has the form: 



m a -)N, = J run 



u=l k=l n v 



k=l m - 



Inserting ( |2.3| ) into ( |2.11| ) we find 

R((cx-)n, (<%+)n) = exp{z E y da;da?'[0+(a:)D 



AT 

'a; — a;'; a+ v )<p-ix' ) — 



-0_(x)D_ + (x - x'; a_,„)^ + (x , )]}^(0 + )^*(0-), (2.12) 
where 0_ is defined on the complex conjugate contour C_ :t~^t — ie and 

D + _(x - x'; a) = -ij dto(q)e iq(x - x ' ) e- iaq , (2.13) 

D_ + (a; -x';a)=i J du{q)e- iq( ~ x -< ] e~ iaq (2.14) 
are the positive and negative frequency correlation functions correspondingly. 



We must integrate over sets (Q) n and (P) n if the distribution of fluxes momenta over 
cells is not fixed. In result, 

r(P) = J D 4 a-(P)d^a + {P)R((a^ N , (a + ) N ), (2.15) 

where the differential measure 

N J& 

DMP) = U^K(P,(*)n) (2.16) 

takes into account the energy-momentum conservation laws: 

N n 



K(P,(a) N ) = f Hd 4 Q u e^ a » Q »5 (4 \P -£Q„). (2.17) 
i 1 

The explicit integration gives that 

K(P,(a) N ) ~ ]j6 (3 \a-a u ) } (2.18) 



u=l u=l 



N 



u=l 



where a is the center of mass (CM) 3- vector. 

To simplify the consideration let us choose the CM frame and put a = (— i/3,0). In 
result, 

roc N N 

K(E, (P) N ) = / J] dE w e^^5(E - £ E v ) (2.19) 
Correspondingly, in the CM frame, 

r(E) = [ Df3 + (E)D[3_(E)R((I3 + ) N , (/3_) JV ), (2.20) 



where 

N d3 

Df3(E) = ft£:K(E,((3) N ) (2.21) 



and R(((3)n) was defined in ( |2.12| ) with otk, v = (— i/3^,0), Ref3k,u > 0, k = +, 



We will calculate integrals over j3 k using the stationary phase method. The equations 
for mostly probable values of (5 k : 

1 <} -K(E,((3 k ) N )= r-| R{(P) N ), k = +,- (2.22) 



K(E,(p k ) N )dp k> „ y R(((3i) N )d(3 k „ 

always has the unique positive solutions ^^(E). We propose that the fluctuations of /3 k 
near /3 k are small, i.e. are Gaussian. This is the basis of the local-equilibrium hypothesis 
||13|| . In this case is the temperature in the initial state in the measurement cell v 

and l/j3 +iV is the temperature of the final state in the v-th measurement cell. 



The last formulation ( |2.15| ) imply that the 4-momenta (Q)n and {P)n can not be 
measured. It is possible to consider another formulation also. For instance, we can 
suppose that the initial set (Q)n is fixed (measured) but (P)n is not. In this case we will 
have mixed experiment: f3- )V is defined by the equation: 



and f3 +tV is defined by second equation in ( |2.22| ). 



Considering limit N — > oo the dimension of cells tends to zero. In this case we are 
forced by quantum uncertainty principle to propose that the 4-momenta sets (Q) and (P) 
are not fixed. This formulation becomes pure thermodynamical: we must assume that 
(/?_) and (/?+) are measurable quantities. For instance, we can fix (/?_) and try to find 
((3+) as the function of total energy E and the functional of (/?_). In this case eqs. ( |2.22j ) 
become the functional equations. 

In the considered microcanonical description the finiteness of temperature does not 
touch the quantization mechanism. Really, one can see from (|2.12|) that all thermody- 
namical information is confined in the operator exponent 

e N(4>i0j) _ j~j j~j g i/ faDi^j (2.24) 

the expansion of which describes the environment, and the "mechanical" perturbations 
are described by the amplitude Z(<p). This factorization was achieved by introduction of 
auxiliary field <p and is independent from the choice of boundary conditions, i.e. from the 
choice of the considered systems environment. 



3 The distribution functions 

In the previous section the generating functional R((P)n) was calculated by means of 
dividing the "measuring device" (calorimeter) on the N cells. It was assumed that the 
dimension of device cells tends to zero (N — > oo). Now we will specify the cells coordinates 
using the Wigner's description |3[] . 

Let us introduce the distribution function F n which defines the probability to find n 
particles with definite momentum and with arbitrary coordinates. This probabilities 
(cross section) are usually measured in particle physics. The corresponding Fourier- 
transformed generating functional can be deduced from ( |2.12 ): 

N , 

F( Z , (/3 + ) N , (p_) N ) = n n^^ 8 x 

u=l v£j 

xZ(0 + )Z*(0_). (3.1) 

The variation of F over z\Aq) generates corresponding distribution functions. One can 
interpret z"Aq) as the local activity: the logarithm of z"Aq) is conjugate to the particles 



number in the cell v with momentum q for the initial {ij = H — ) or final (ij = — h) states. 
Note that z^(q)(f>*(q)(pj(q) can be considered as the operator of activity. 

The Boltzman factor e - ^^^ can be interpreted as the probability to find a particle 
with the energy e(q) in the final state (i = +) and in the initial state (i = —). The total 
probability, i.e. the process of creation and further absorption of n particles, is defined 
by multiplication of this factors. 

The generating functional ([O]) is normalized as follows: 

F(z = l,(P))=R((J3)), (3.2) 

F{z = 0, 03)) = \Z(0)\ 2 = M<f>±)U±=o (3.3) 

Where 

i? o (0 ± ) = Z(0 + )Z*(0_) (3.4) 

is the "probability" of the vacuum into vacuum transition in presence of auxiliary fields 
4>±. The one-particle distribution function 

F 1 {((3 + ) N , (/?_)„; q) = jJ-^FU=o = 

= {^(g)e^ e ^/ 2 }{^(g)e- /3 ^ ( " )/2 } J Ro(0 ± ) (3.5) 

describes the probability to find one particle in the vacuum. 
Using definition (|2.4j) , 

F 1 ((/3 + ) w ,(/5_) w ;g) = J dxda;V 9(a; -^ ) e-' 3 - e(9) }0 J (x)^(x') J Ro(0±) = 

= J dYidye^e-^^UY + y/2)UY - y/2)R (<P ± )}. (3.6) 

We introduce using this definition the one-particle Wigner function W\ @: 

F^P+^iP^N-q) == J dYW 1 ((f3 + ) N ,((3-) N -Y,q). (3.7) 

So, 

W 1 ((P + ) N ,(P^) N ;Y,q) = J dye^e-^^UY + y/2)M Y -V/ 2 )M4>±)- (3.8) 

This distribution function describes the probability to find in the vacuum particle with 
momentum q at the point Y in the cell v 

Since the choice of the device coordinates is in our hands it is natural to adjust the 
cell coordinate to the coordinate of measurement Y: 

WxiiP+UM-.WXq) = J dye^e-^ Y ^}4> t (Y + y/2)4j(Y - y/2)Ro(<j>±). (3.9) 



This choice of the device coordinates lead to the following generating functional: 

F(z,(3) =exp{ i J dydY[4> + (Y + y/2)D + _(y;(3 + (Y),z)^(Y-y/2)- 

+ y/2)D_ + (y; /3_(Y), z)MY - y/2)]}R ((f> ± ), (3.10) 

where 



D + 4y;P + (Y),z) = -i J du(q)z + 4Y,q)e^e-^ Y ^\ (3.11) 
D„ + (y;(3 + (Y),z) = i f du(q)z^ + (Y,q)e-^e-^ Y ^ (3.12) 



are the modified positive and negative correlation functions ( |2.13j ), ( |2.14| ). 



The inclusive, partial, distribution functions are familiar in the particle physics. This 
functions describe the distributions in presence of arbitrary number of other particles. 
For instance, one-particle partial distribution function 

-f3i(Y)e(q) . 

= {27r) 3 e{q) J dye^^Y + y/2)^-(y - y/2)R{</> ± , (/3)), (3.13) 

where eq.(|3.2|) was used. 

The mean multiplicity (Y, q) of particles in the infinitesimal cell Y with momentum 

q is 

nij (Y,q) = J dq lnF(z,(/3))U. (3.14) 

If the interactions among fields are switched out we can find that (omitting indexes): 

n ( Y > = e p(Y) qo _x ' * = e (?) > °- ( 3 - 15 ) 
This is the mean multiplicity of black-body radiation. 



4 The closed-path boundary condition 

The developed in Sec. 2 formalism allows to introduce the more general boundary condi- 
tions instead of fl2.5|) . Considering the probability R which has the double path integral 
representation we will introduce integration over closed path. This allows to introduce 
the equality: 

I da^+d^ + - $_d M $_) = 0, (4.1) 

•'Too 

as the boundary condition, where o"oo is the infinitely far hypersurface. The general 
solution of this equation is: 

$ ± ( C T 00 ) = $(a 00 ) (4.2) 



where $(<7oo) is the "turning-point" field. The result of this changing of boundary condi- 
tion was analyzed in [12| for the case of uniform temperature distribution. 

In terms of .S-matrix the field ^((Too) represent the background flow of mass-shell 
particles. We will propose that the probability to find a particle of the background flow 
is determined by the energy-momentum conservation law only. In another words, we will 
propose that the system under consideration is surrounded by the black-body radiation. 

Presence of additional flow will reorganize the differential operator exp{iV(0j0,,)} only 
and new generating functional R cp has the form: 

R cp (a + ,a_) = e N ^R (<P ± ). (4.3) 



The calculation of operator iV(0j0,,) is strictly the same as in |12| . Introducing the cells 
in the Y space we will find that 



Nifah) = J dYdyUY + y/2)nij(Y, y)<P 3 (Y - y/2), (4.4) 
where the occupation number carries the cells index Y: 

n l3 (Y, y) = J du^e^n^Y, q) (4.5) 

and (g = e{q)) 

n ++ (Y,q ) = n__(F,g ) = n(Y, + /?_)|g b |/2) = e(/3++/3 _ )( y )ko | /2 _ p (4-6) 

n + _(Y,q ) = Q(q )(l + h(Y,(3 + q )) + Q(-q )n(Y,-P„q ), (4.7) 

n-+(y,q )=n + -(y,-q ). (4.8) 

For simplicity the CM system was used. 

Calculating R perturbatively we will find that 

R cp {(3) = exp{-iV(-ij+)+iV{-ij.)} x 

x exp{* J dYdy\ji{Y + y/2)G tJ (y, (f3(Y))jj(Y - y/2)} (4.9) 

where, using the matrix notations, 





and 



iG(q, ((3(Y))) = «^+- j + 

\ q 2 —m?—ie I 

+to6tf-m*)( n ^ 1±± f m ^ n((3 + (Y)\ qo \)a + (M\ Uw) 

+ 2,5(q m ) ^ n(/3 _ (y)|go|)a _ (/3 _ ) n{ iJ±±Mm M J' ^ 10 ) 

g ± {j3) = _ e /5(l<?0|±<?0)/2_ ( 4U ) 

Formally this Green functions obey the standard equations in the y space: 

(&-m 2 ) y G u = 6(y), 

(<9 2 -m 2 ),G^ = 0, i^j (4.12) 

since ^CToo) ^ reflects the mass-shell particles. But the boundary conditions for this 
equations are not evident. 



5 Concluding remarks 



One can not expect the evident connection between the above considered and Zubarev's 
jn| approaches. The reason is as follows. 

In Zubarev's theory the "local-equilibrium" hypothesis was adopted as the boundary 
condition.lt is assumed that in the suitably defined cells of a system at a given temper- 
ature distribution T(x,t) = l/(3(x,t) where (x,t) is the index of the cell, the entropy is 
maximum. The corresponding nonequilibrium statistical operator 

Rz ^ e ~ J d 3 x?T 00 ( 51 ) 

describes evolution of a system. Here T^ u is the energy-momentum tensor. It is assumed 
that the system "follows" to (3(x,t) evolution and the local temperature T(x,t)is defined 
as the external parameter which is the regulator of systems dynamics. For this purpose 



the special ie-prescription was introduced fT3|j . 



The KMS periodic boundary condition || [TO]] can not be applied HIT , [12| and by this 
reason the decomposition: 

P&t)=p + p 1 (2,t) (5.2) 
was offered in the paper |J. Here (3q is the constant and the inequality 

Po» \Pi(£,t)\ (5.3) 

is assumed. Then, 

R z ~ e -MH +v + B) (5 _ 4) 

where Hq is the free part of the Hamiltonian, V describes the interactions and the linear 
over Pi/fio term B is connected with the deviation of temperature from the "equilibrium" 
value l//?o- Considering V and B as the perturbations one can calculate the observables 
averaging over equilibrium states, i.e. adopting the KMS boundary condition. Using 
standard terminology || one can consider V as the "mechanical" and T as the "thermal" 
perturbations. 

The quantization problem of operator ( |5.4fj is connected with definition of the space- 
time sequence of mechanical (V) and thermal (B) excitations. It is necessary since the 
mechanical excitations give the influence on the thermal ones and vice versa. It was 
assumed in || that V and B are commuting operators, i.e. the sequence of V- and 
.B-perturbations is not sufficient. This solution leads to the particles propagators renor- 
malization by the interactions with the external field /3(x,t) even without interactions 
among fundamental fields. (Note absence of this renormalizations in our formalism.) 

In the operators V and B are noncommuting ones and -B-perturbations were 
switched on after V^-perturbations. In this formulation the nondynamical renormaliza- 
tion are also present but it is not unlikely that they are canceled at the very end of 
calculations |15| . 



This formulation with (3(x, t) as the external field remained the old, firstly quantized, 
field theory in which matter is quantized but fields are not. It is known that consistent 



quantum field theory requires the second quantization. Following to this analogy, if we 
want to take into account consistently the reciprocal influence of V- and -B-perturbations 
the field j3(x, t) must be fundamental, i.e. must be quantized (and the assumption of paper 
|| becomes true). But it is evidently the wrong idea in the canonical Gibbs formalism. So, 
as in the firstly quantized theory, the theory with operator ( |5.1| ) must have the restricted 
range of validity [[L3[ . 



Therefore, we must reduce our formalism just to the hydrodynamical accuracy to find 
the connection with Zubarev's approach. There is the another side of this question. The 
offered formalism is able to describe an arbitrary nonequilibrium process since it based 
on the ^-matrix, i.e. on the strict field-theoretical description. But the mechanism of 
irreversibility is not clearly seen: the generating functional Ro{4>±) is described by the 
closed-path motion in the functional space, i.e. formally is the time-reversible quantity. 
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